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We consider the effect of optical depth of the 2 3 S level on the nebular recom- 
bination spectrum of He I for a spherically symmetric nebula with no systematic 
Q ■ velocity gradients. These calculations, using many improvements in atomic data, 

can be used in place of the earlier calculations of Robbins. We give representative 
Case B line fluxes for UV, optical, and IR emission lines over a range of physical 
conditions: T = 5000 - 20000 K, n e = 1 - 10 8 cm" 3 , and r 3889 = - 100. A 
^ ' FORTRAN program for calculating emissivities for all lines arising from quantum 

■ levels with n < 10 is also available from the authors. 

We present a special set of fitting formulae for the physical conditions relevant 
to low metallicity extragalactic H II regions: T = 12, 000 — 20, 000 K, n e = 
1 — 300 cm -3 , and r 3889 < 2.0). For this range of physical conditions, the 
Case B line fluxes of the bright optical lines 4471 A, 5876 A, and 6678 A, are 
changed less than 1%, in agreement with previous studies. However, the 7065 
A corrections are much smaller than those calculated by Izotov & Thuan based 
on the earlier calculations by Robbins. This means that the 7065 A line is a 
better density diagnostic than previously thought. Two corrections to the fitting 
functions calculated in our previous work are also given. 



Subject headings: atomic data — ISM:general — ISM: abundances — cosmology: 
observations — cosmological parameters 
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1. Introduction 

In a previous paper (Benjamin, Skillman, & Smits 1998, Paper I), we presented a calcu- 
lation of the nebular recombination line spectrum of He I which combined the atomic data 
and recombination cascade of Smits (1991; 1996) with the collision strengths data of Sawey 
& Berrington (1993). The resulting Case B line fluxes were tabulated, fitting functions were 
provided, and estimates of uncertainties in atomic data on the helium abundance estimates 
were also given. 

An additional uncertainty in estimating helium abundances is the effect of radiative 
transfer. By assuming Case B, we have assumed that most nebulae have a large optical depth 
in transitions arising from the 1 1 S level. However, a significant population of electrons may 
also build up in the metastable 2 3 5 l level. Absorption and reprocessing of photons from 
the 2 3 S level to upper levels will alter the recombination cascade and resulting helium 
abundance estimates. 

The last full treatment of the recombination cascade combined with optical depth was 
provided by Robbins (1968), and this work remains the classic reference for most researchers. 
Additional work on this problem has been done by Almog & Netzer (1989), Proga, Mikola- 
jewska, & Kenyon (1994), and Sasselov & Goldwirth (1995). However, these papers, while 
they considered the radiative transfer problem, used models of the helium atomic physics 
with an accuracy less than the 1% level necessary for meaningful cosmological constraints. 
As a result, while they were useful in demonstrating the relative effect of radiative transfer; 
the absolute effect has not been addressed since the improvements in the atomic data. We 
address this problem here. 

This work is of particular interest to the derivation of the primordial helium abundance. 
Accuracy in the atomic physics of better than 1% is required to produce abundance estimates 
that provide meaningful constraints on models of Big Bang nucleosynthesis. Currently, ob- 
servational studies are claiming accuracies on order of less than 1% (e.g., Y p = 0.2452 ± 
0.0015, Izotov et al. 1999; Y p = 0.2345 ± 0.0026, Peimbert et al. 2000). In Paper I we per- 
formed a Monte Carlo exercise to estimate the uncertainty in the theoretical emissivities due 
to uncertainties in the input atomic data and obtained an estimate of 1.3 to 1.5%. Recently, 
Olive & Skillman (2001) have demonstrated that minimization routines used to solve for 
physical conditions and He abundances simultaneously can result in erroneous solutions due 
to the degeneracies of the dependencies of the helium emission lines on the various physical 
conditions (i.e., temperature, density, optical depth, and underlying absorption). Obviously, 
using accurate atomic physics in these calculations reduces the systematic uncertainties in 
the abundance calculations. 
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2. Calculations 



For this work, we consider the radiative transfer of a spherically symmetric nebula with 
no systematic expansion. We assume that the recombination is in the case B limit, so that 
the only radiative transfer effects are those that arise from the metastable 2 3 S level. The 
assumption of no systematic expansion is justified by detailed studies of the velocity fields 
of giant extragalactic HII regions (e.g., Yang et al. 1996, and references therein) which show 
them to be dominated by turbulent, and not organized, motions. We parameterize our cal- 
culations in terms of the line center optical depth of the 3889 A line, T3889 = n(2 3 S )asssdRs, 
where n(2 3 S ) is the density of atoms in the metastable state. 1 The cross section at line 
center is given by 



a = (6.73 x 10" 31 cm" 2 ) A\ 



1-1/2 



(1) 



where A(s _1 ) is the spontaneous transition probability; \^ is the wavelength in Angstroms; 
the gas temperature is T 4 = T/(10 4 K); m He = 4.0, and the "turbulent" velocity, £, is 
typically 40 km s" 1 . 

For a nebula of Stromgren radius, R s = (30 V c )Qis n H^-> the optical depth is 



r = 9.68 x 10- 14 Q^ 3 n^V 7 AT 4 - 1/2 A^ , (2) 

where Q 48 = Qh/(10 48 cm" 3 ) is the luminosity of hydrogen ionizing photons, n H (cm" 3 ) is 
the hydrogen particle density, A 6 = A/(10 6 s" 1 ) is the spontaneous transition probability, 
and p_ 7 = (n(2 3 S )/n He+ )/lO~ 7 is the population in the metastable 2 3 S level. A helium 
abundance of 10% by number is assumed. Turbulent velocities will lower this optical depth, 
and can be accounted for by replacing T 4 with T 4 + m/f e ^ 2 /166, where £ is in km s -1 . For low 
densities, the population in the 2 3 S level is in the range of p_7 = 1-4 (See Table 8 in Paper 
I.) For fixed T and p, the variation of optical depth with ionizing luminosity and particle 
density is shown in Figure 1. 

For a spherically symmetric nebula, the mean probability of escape of a photon with a 
frequency characterized by x = (u — v )/Au is given by (Cox & Mathews 1969) 



1 Note that in the text we refer to emission lines using their wavelength in air. However, the tables use 
the vacuum wavelengths derived from the energy levels described in Paper I. 
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Fig. 1. — Dependence of optical depth in 3889 A line center as a function of ionizing lu- 
minosity of source, Qh, and hydrogen particle density, n#. This is assuming a fixed he- 
lium abundance of 10% (by number), a temperature of 10 4 K, no turbulent velocity, and 
n(2 3 S )/riHe+ = 10~ 7 . Using a typical turbulent velocity of £ = 40 km s -1 would reduce the 
optical depth by a factor of 6.2. The scaling of optical depth with these quantities is given 
in equation 1 and 2. 



-2t x 



(3) 



where v is the line center frequency and Au = {v$/c) ^j2kT /m He + £ 2 . 

Integrating this escape probability over the entire line profile yields a total escape prob- 
ability of e(r) = ^= J^° oo p(T x )e~ x2 dx. As noted in Osterbrock (1989), this can be approxi- 
mated to within 10% with the function e(r) = 1.72/(1.72 + r). To get a more exact result, 
we have evaluated this integral numerically, with a precision of better than 0.01%. 

The calculation of statistical equilibrium given in equation (4) of Paper I is then simply 
modified by replacing the spontaneous decay rates A with e(r)A, yielding a revised set of 
level populations and emissivities. In the calculations here, we used a model atom with 
individual levels up to n = 20; otherwise the calculations are the same as in Paper I. The 
relative optical depths of lines in the n 3 P-2 3 S series are given in Table I. 
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3. Results 

3.1. Tables, Fitting Formulae and Comparisons 

We have calculated the emissivity, j = j (n, T, r 3889 ), of recombination lines for the same 
temperature and density grid as in Paper I, but for optical depths up to r 38 89 = 100. The 
ratio of the line emissivity for optical depth r 388 9 to the emissivity for zero optical depth 
is given by optical depth factor /We (73889) = j( n ,T,r 38 sg)/ j(n,T,0). Figure 2 shows how 
the optical depth factor of several lines varies as a function of optical depth for a case with 
n e = 10 2 cm -3 and T = 10000 K and 20000 K. Although we include the collisional coupling 
between the triplet and singlet levels, only the triplet lines show noticeable differences. Even 
for densities as high as n e = 10 8 cm -3 , the effect of optical depth in the 2 3 S level on the 
singlet lines is small, less than 0.4% in the worst case. 




20 40 60 80 100 



T 3889 

Fig. 2. — Optical depth correction factor for selected optical lines for cases with a particle 
density of tin = 100 cm -3 and two different temperatures. 

Table 2 shows all optical and IR lines that change by more than 5% over the range 
r 3889 = — 100 for the case with n e = 10 2 cm -3 and T = 10000 K. This provides a guide to 
which lines are most affected by optical depth. From Figure 2, note the strong dependence of 
the 7065 A line on both r 3889 and temperature. The decrease in /7065 (73889) with increasing 
temperature (at constant density) arises due to a decreasing population in the 2 3 S level 
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caused by the increasing collisional excitations out of that level. We also provide polynomial 
fit coefficients in Table 3 for these lines that approximate the calculated results to within 
10%, and better than 1% for a few selected optical lines at a single temperature and density. 

Unfortunately, the optical depth factor is sometimes a complex function of density, 
temperature, and r 3889 , so no simple functional form was found that could provide a good 
(i.e., better than 10%) fit for fu n e{ T 3&m) over our full three dimensional parameter space. 
This difficulty is illustrated by Figure 3 which shows /7065 at t 388 9 = 100. As a result, we 
have written a program which does a linear interpolation on our grid of results to estimate 
the emissivities. This program may be obtained from the authors. 2 



Countours of f 7 065( T 3S89= 1 00) 




2 4 6 8 

log, n (cm -3 ) 



Fig. 3. — Value of optical depth correction factor for the 7065 A line at T3 889 = 100 for a 
range of temperatures and densities. There is a large range of variation that is not well fit 
by a simple function. As a result, those desiring accurate estimates of these lines should use 
the FORTRAN program available from the authors. 

In Figure 4, we have compared the results of our new calculations with the results of 
Robbins (1968) for a density of n e = 100 cm -3 . Robbins' results are shown as filled symbols, 
ours are shown as open symbols, and a solid line shows the result of the fitting formulae 
in Table 3. This figure shows differences between the relative line strengths for different 



2 Program and input data file may be obtained from the website http://wisp.physics.wisc.edu/~benjamin. 
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temperatures, particularly for A7065. This demonstrates again that the radiative transfer 
effects are dependent on other factors beside t 3889 . As a result, the fitting formulae given in 
Table 3 should be used with caution, since they are only for a single choice of temperature 
and density. 

p 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l 

10,000 K(R68) 
10,000 K(BSS) 
20,000 K (R68) 
20,000 K (BSS) 
BSS fit — 
IT(98) lit 
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Fig. 4. — Comparison of recent work to earlier work of Robbins (1968) and the Izotov & 
Thuan (1998) fits to the Robbins (1968) calculations. The dotted vertical lines indicate the 
range of values of r(3889) that are important for the giant extragalactic HII regions that 
are used for determining the primordial helium abundance. Note that the range in intensity 
ratio on the y-axis changes signficantly from panel to panel. In some cases (e.g., A3889) the 
differences in the results for the two different temperatures are small compared to the y-axis 
range, and thus, the circular and square symbols are over-plotted. 

Second, although Robbins only reported results on a very coarse grid, it is instructive 
to note the difference in / 38 89 as a function of r 388 9. Our calculated A3889 line strengths 
do not decrease as strongly with increasing T3 88 g. Since A3889 is essentially a resonance 
line, most A3889 absorptions are followed directly by a A3889 emission, with the result that 
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the line strength falls off relatively slowly with increasing T3889 (less than 50% at r 38 89 = 
fO). The absolute value of ^3889 as a function of r 38 89, therefore, depends on the formalism 
used in the calculation. Robbins (1968) used an "exact redistribution" function from Unno 
(1952), while here we use complete redistribution. In the complete redistribution case, 
the absorbed wavelength and re-emitted wavelength are completely unrelated. In addition, 
Robbins (1968) used a numerical integration of the equations of radiative transfer while our 
result is (almost) analytical, and therefore has much higher numerical accuracy. Thus, we 
believe our calculations should be the more accurate. 

A significant advantage to using an escape probability formalism as opposed to calcu- 
lating the "exact" solution as in Robbins (1968) is that it allows one to quickly change the 
characteristics of the radiative transfer. Since these effects depend upon the exact geometry 
of the system, which will always be poorly known, our ability to correct line fluxes for modifi- 
cations due to self-absorption will always be limited. In particular, density inhomogeneities, 
ionization inhomogeneities, large scale velocity structure, and planar versus spherical geom- 
etry can all effect the resulting radiative transfer to differing degrees. Our calculations also 
do not include the effects of dust absorption or the effects of the hydrogen emission line at 
3889.1 A. 

In lieu of calculating a large number of grids of differing geometries and density /velocity 
structure, we have attempted to estimate the robustness of our results using two alternate 
forms of the escape probability corresponding to differing geometries and scattering as- 
sumptions. For the case with T = 10 4 K and n e = 10~ 2 cm" 3 , we have recalculated the 
recombination cascades using the simple form e = (1 — exp(— r))/r for a plane parallel ge- 
ometry (Rybicki 1984; Almog & Netzer 1989); the second form comes from integrating the 
local escape probability given by equation 1.2 in Avrett & Hummer (1965). We find that 
for the restricted range of optical depth of interest for primordial helium abundance stud- 
ies, the difference in all the bright optical lines is less than 1%. For the extended range to 
73889 — 100, all the triplet lines that are not part of the 2 3 S' — n 3 P series are modified by 7% 
at the most. Lines in the 2 3 S — n 3 P series are affected much more. Like other lines, they 
are modified indirectly by the change in the 2 3 S' level population (which is a small fractional 
change in the line flux) , but also directly by the reprocessing of photons in this series. 

Another important effect for radiative transfer is the presence or absence of velocity 
gradients. Our treatment implicitly assumes that the expansion rate of a nebula is signif- 
icantly less than the intrinsic Doppler width of the line, i.e., w = v exp jv^ ovv << 1, where 
v dopp — ^/lQQTi/rriHe + C 2 as above. For extragalactic HII regions, we are generally justified 
in assuming that w = 0. However, using the results of Robbins (1968), one can estimate 
the effects of expansion velocity for an uniform density sphere. We find that the effects of 
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expansion can be approximated by replacing r 388 9 with T 388 g ie // — T3889/(l + 0.4iu). This 
approximation is good to within 5% up to r 388 9 = 30 and 20% for r 3889 = 70. 

3.2. Relevance to Primordial Helium Determinations 

Part of the motivation of this work was to estimate accurate radiative transfer correc- 
tions in determinations of helium abundances in low metallicity extragalactic HII regions. 
Thus, as in Paper I, we have calculated optical depth correction factors over a limited grid 
of physical parameters: T = 12, 000 - 20, 000 K, n e = 1 - 300 cm" 3 , and r 3889 < 2.0. The 
upper bound on the optical depth is chosen for two reasons. First, most observations of ex- 
tragalactic HII regions find optical depths less than 2.0 (e.g., Izotov & Thuan 1998). This is 
due to the large (> 25 km s _1 ) velocity dispersion typical of these objects (e.g., Melnick et al. 
1988, and references therein). Second, if an object requires a large optical depth correction to 
derive the helium abundance, it is probably an inherently uncertain abundance calculation, 
and may not therefore be useful in determining the primordial helium abundance. 

The upper bound on density follows from the same principle, i.e., as the correction 
for collisional excitation becomes large, the uncertainty in that correction must also grow, 
rendering the abundance determination less certain. While the density range covers the 
entire range of densities considered for NGC 346 by Peimbert et al. (2000) and almost all 
of the targets studied by Izotov & Thuan (1998), Peimbert, Peimbert, & Luridiana (2001) 
have derived much higher densities for several of the targets in the Izotov & Thuan (1998) 
sample. If it proves necessary to consider a grid with a larger density range, it is possible to 
produce one with the provided programs. 

We find that the optical depth correction factors for the bright optical lines 4471 A and 
5876 A are less than 1% in this limited range. It is the corrections to 7065 A and 3889 A 
that are of most concern. In Table 4, we give fitting functions that allow one to estimate the 
helium abundance given an observed line ratio of r/ ine = Iune/^H/3- The helium abundance 
by number for a given line is given by 

Vline ^Zine"7 7 7 

Jline{n e , t, T 3889 

The function Fi ine (n e) t) = At Bo+Bine is taken from Paper I, while the optical depth function, 
fiine(n e ,t,T 38 s9) = 1 + (r / 2)[a + (b + bin e + b2nl)t] i s calculated here for the limited ranger < 
2.0. For the 7065 A line however, it is necessary to use a different functional form. The Hf3 
emissivity is fit by tHp{n e) t) = 4:7ijHf3/n e nH+ = 12.450 x 10~ 26 t~ ' 917 with an maximum error 
of 0.3% over our limited parameter space. The 7065 emissivity at r = is fit by 67065 {n e , t) = 



(4) 
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^j 7om /n e n He+ =C + (D + Dm e + D 2 n 2 e )t, where (C, D , D u D 2 ) = (3.4940, -0.793, 1.50 x 
10~ 3 , — 6.96 x 10~ 7 ) x 10 -26 , and the optical depth function /Vo65(^ e , t, r 3889 ) is the same 
as above with parameters (a,6 AA) = (0.359,-3.46 x 10~ 2 , -1.84 x 10~ 4 , 3.039 x 10~ 7 ). 
Then the helium abundance by number using the 7065 line is given by 



e H/3 (n e ,t) 

1/7065 — ?>065 7 777 7 7 7 \° ) 

^7065 {n e , t) /7065l^e, £, ^3889 ) 

These functions fit 2/7065 within 2.2% over our limited range. 

Figure 5 shows the same functions as in Figure 4 but only for the low values of r 3889 
relevant to estimating the primordial helium abundance. We include the fits to the optical 
depth correction factors contained in Table 4, and refer to them as "BSS limited fit". A 
striking feature in Figure 5 is the large deviation of the Izotov & Thuan fitting formula from 
our newly calculated values for the 7065 A line. Since the Izotov & Thuan fit is based on 
the very coarse grid of values given by Robbins, it does not accurately represent the results 
at low T3889 where it is actually being applied. The parabolic form of the fitting function 
is not appropriate to the physics, as / increases nearly linearly with r 3889 for low values of 
r 3889 , and turns over only at relatively high values of r 3889 . 

Evidence of this can be seen by considering the predictions of the Izotov & Thuan 
functions at r 3889 = 0.5. Figure 5 shows that at r 3889 = 0.5 the 3889 A line has decreased by 
roughly 3% while the Izotov & Thuan fit for the 7065 A line has increased by roughly 25%. 
At 20,000 K, the ratio of the emissivities for the A3889 and A7065 lines is 3.85, which, after 
correcting for energy per photon, means that there are roughly 2.1 A3889 photons emitted 
for each A7065 photon. Since the increase in the 7065 A line with increasing r 3889 comes 
from the conversion of A3889 photons, the 7065 A line should not increase faster than 2.1 
x the 3889 A line decreases. Thus, in the range of interest, the Izotov & Thuan fit has 
overestimated the radiative transfer effects on the 7065 A line by roughly a factor of four at 
r 3889 = 0.5, and higher at lower values of r 3889 . 

This difference has important consequences. Olive & Skillman (2001) explored various 
possibilities for systematic errors that could result by using minimization routines to solve 
for physical conditions and He abundances simultaneously. Since the 7065 A line plays the 
main role as the density diagnostic in these minimization routines, overestimating the 7065 
A dependence on r 3889 leads to errors in the density estimates (which directly affect the 
helium abundance determinations). The good news is that the 7065 A line dependence on 
r 3889 is much less than previously thought, so that using the 7065 A as a density diagnostic 
will increase in diagnostic power. That is, the degeneracies between the 3889 A and 7065 
A lines will be decreased. This is important since the 3889 A line is subject to systematic 
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Fig. 5. — A plot similar to that shown in Figure 4, zooming in on the range of r(3889) from 
to 3. Note that the y-axis ranges are different from those of Figure 4, and, again, change 
from panel to panel. The striking feature in this figure is the large difference between the 
new calculated values for A7065 and the interpolation fit to the Robbins (1968) calculations 
derived by Izotov & Thuan (1998). For low values of r(3889), the new calculations show 
that the amplification of the A7065 line is much lower than would have been inferred from 
the Izotov & Thuan (1998) fit. This implies that the A7065 line is a better density diagnostic 
than previously thought. As in Figure 4, the differences in the results for the two different 
temperatures for A3889 are small compared to the y-axis range, and thus, the circular and 
square symbols are over-plotted. 

observational uncertainties, most notably separation from the 3889 HI line. 

The error estimation work of Olive & Skillman (2001) will need to be revisited with 
the newly calculated fitting formulae. We recommend that all future helium abundance 
calculations from minimization routines take advantage of the new fitting formulae presented 
in Table 4. The consequence of our results on the determination of the primordial helium 
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abundance is not yet clear. Only about one third of the targets in the sample of Izotov 
& Thuan (1998) showed evidence for detectable radiative transfer effects. Moreover, since 
the main effect comes through the densities calculated from the 7065 A line, the effect is 
likely to be small in most cases. The total effect will need to be calculated on a case-by-case 
basis. However, the sense of the correction will be that higher 7065 A line strengths will 
be interpreted as being due to higher densities, not higher values of r 38 89- As can be seen 
using the helium abundance fitting functions in Table 9 of Paper I, this will result in a lower 
helium abundance for objects for which radiative transfer corrections have been applied. 

We would like to thank John Mathis for comments and encouragement, and Manuel 
Peimbert and Stuart Pottasch for pointing out errors in Paper I. We would also like to 
acknowledge the support of NASA Astrophysics Theory Grant NAG5-8417 (RAB) and NASA 
LTSARP grant NAG5-9221 and the University of Minnesota (EDS). 

A. Corrections to Paper I 

Two corrections to the fitting formulae provided in Paper I should be noted. First, 
the fitting formula for primordial helium abundance in Table 9 for the 6678 line should be 
given by / = 2.58t a249 - ne2 - 0xKr4 . Second, the values of fluxes in Tables 5-7 for Brackett 7 
emissivities were in error. The values of the emissivities in units of 10~ 27 ergs s _1 cm 3 for 
T = (5000,10000,20000 K) should be (7.25,3.44,1.54), (7.20,3.41,1.53), and (6.99,3.33,1.51) 
for n e = 10 2 , 10 4 ,and 10 6 cm -3 . The corresponding fitting functions (good to within 8%) 
are 3.39 x lO" 27 *" 1 - 118 , 3.37 x lO" 27 *" 1 - 116 , and 3.29 x lO" 27 *" 1 - 104 . 
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Table 1. Optical depth ratios of 2 3 S-n 3 P series' 



n u 


A(A) 


T 


n u 


A(A) 


T 


2 


10833 


2328.35 


12 


2654 


0.86 


3 


3890 


100.00 


13 


2646 


0.66 


4 


3189 


32.77 


14 


2639 


0.51 


5 


2946 


14.68 


15 


2634 


0.42 


6 


2830 


7.88 


16 


2630 


0.34 


7 


2765 


4.74 


17 


2627 


0.30 


8 


2724 


3.08 


18 


2624 


0.23 


9 


2697 


2.12 


19 


2621 


0.21 


10 


2678 


1.53 


20 


2619 


0.18 


11 


2664 


1.10 









a Normalized to r 389 o = 100. 
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Table 2. Optical Depth Correction //m e (r 3889 ) for He I Emission Line Intensities a 



A (A) 


n u 


/ 


n 9 


n — in 


T — DU 


- — i nn 

T — 1UU 


2946 


5 


1 


0.96 


0.81 


0.40 


0.29 


3189 


4 


1 


0.93 


0.71 


0.29 


0.19 


3733 


7 





1.00 


1.01 


1.04 


1.06 


3868 


6 





1.00 


1.02 


1.09 


1.12 


3890 


3 


1 


0.90 


0.64 


0.22 


0.14 


4122 


5 





1.01 


1.06 


1.20 


1.25 


4472 


4 


2 


1.00 


1.01 


1.05 


1.06 


4714 


4 





1.04 


1.19 


1.54 


1.64 


5877 


3 


2 


1.01 


1.03 


1.07 


1.08 


7067 


3 





1.30 


2.10 


3.52 


3.81 


7162 


10 


1 


1.00 


1.04 


1.24 


1.37 


7300 


9 


1 


1.00 


1.06 


1.32 


1.47 


7501 


8 


1 


1.01 


1.09 


1.45 


1.63 


7818 


7 


1 


1.02 


1.14 


1.63 


1.84 


8363 


6 


1 


1.04 


1.24 


1.91 


2.16 


9012 


10 


1 


1.00 


1.04 


1.24 


1.37 


9231 


9 


1 


1.00 


1.06 


1.32 


1.47 


9466 


5 


1 


1.10 


1.46 


2.41 


2.69 


9555 


8 


1 


1.01 


1.09 


1.45 


1.63 


9705 


7 





1.00 


1.01 


1.04 


1.06 


10075 


7 


1 


1.02 


1.14 


1.63 


1.84 


10671 


6 





1.00 


1.02 


1.09 


1.12 


10833 


2 


1 


1.02 


1.07 


1.17 


1.19 


10999 


6 


1 


1.04 


1.24 


1.91 


2.16 


12531 


4 


1 


1.26 


2.03 


3.56 


3.93 


12850 


5 





1.01 


1.06 


1.20 


1.25 


12988 


5 


1 


1.10 


1.46 


2.41 


2.69 


14492 


10 


1 


1.00 


1.04 


1.24 


1.37 


15067 


9 


1 


1.00 


1.06 


1.32 


1.47 


15951 


8 


1 


1.01 


1.09 


1.45 


1.63 


17006 


4 


2 


1.00 


1.01 


1.05 


1.06 


17382 


10 


1 


1.00 


1.04 


1.24 


1.37 


17455 


7 


1 


1.02 


1.14 


1.63 


1.84 


18214 


9 


1 


1.00 


1.06 


1.32 


1.47 


19523 


8 


1 


1.01 


1.09 


1.45 


1.63 


19550 


4 


1 


1.26 


2.03 


3.56 


3.93 


20428 


6 


1 


1.04 


1.24 


1.91 


2.16 


21128 


4 





1.04 


1.19 


1.54 


1.64 



a These correction factors for case with n e = 10 2 cm 3 
and T = 10 4 K. Only lines from n u < 10 with 1.05 > 
/( T 3889) •> 0.95 over the range T3889 = [0, 100] are included. 
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Table 3. Fitting function for /«i„ e (r 3889 ) for n e = 10 2 cm" 3 and T = 10 4 K. a 



A (A) 


fo 


10 2 /i 


10 4 /2 10 6 /3 10 8 / 4 10 10 /5 Max Err (%) 


2946 


1.000 


-2.054 


2.398 -1.066 ■■■ ■■■ 3.2 


3189 


1.000 


-3.479 


7.466 -7.871 3.082 ■■■ 4.8 


3733 


1.005 




5.4 


3868 


1.010 




9.6 


3890 


1.000 


-4.983 


16.086 -28.586 25.038 -8.415 5.0 


4122 


1.003 


0.288 


3.7 


4472 


1.000 


0.125 


-0.067 •■■ ••■ ■■■ 0.3 


4714 


1.003 


1.485 


-0.896 ■■■ ■■■ ■■■ 4.3 


5877 


1.000 


0.192 


-0.119 ■■■ ■■■ ■■■ 0.7 


7067 


1.000 


15.008 


-46.190 80.374 -69.603 -23.231 1.3 


7162 


0.999 


0.385 


1.4 


7300 


1.000 


0.509 


2.3 


7501 


1.002 


0.687 


3.9 


7818 


1.005 


0.947 


6.8 


8363 


1.001 


2.231 


-1.111 ■■■ ■■■ ■■■ 2.7 


9012 


0.999 


0.385 


1.4 


9231 


1.000 


0.509 


2.3 


9466 


1.006 


3.778 


-2.187 ••• ••■ •■■ 6.8 


9555 


1.002 


0.687 


3.9 


9705 


1.005 




5.4 


10075 


1.005 


0.947 


6.8 


10671 


1.010 




9.6 


10833 


1.005 


0.242 


5.2 


10999 


1.001 


2.231 


-1.111 ■■• ••■ •■■ 2.7 


12531 


1.007 


10.565 


-14.802 7.238 ■•■ •■■ 6.0 


12850 


1.003 


0.288 


3.7 


12988 


1.006 


3.778 


-2.187 ■■■ ■■■ ■■■ 6.8 


14492 


0.999 


0.385 


1.4 


15067 


1.000 


0.509 


2.3 


15951 


1.002 


0.687 


3.9 


17006 


1.006 




5.1 


17382 


0.999 


0.385 


1.4 


17455 


1.005 


0.947 


6.8 


18214 


1.000 


0.509 


2.3 


19523 


1.002 


0.687 


3.9 


19550 


1.007 


10.565 


-14.802 7.238 ■■• ■■■ 6.0 


20428 


1.001 


2.231 


-1.111 ■■• ••■ •■■ 2.7 


21128 


1.003 


1.485 


-0.896 ■■• ••■ •■■ 4.3 


a /(T3889) = fty 


f/ir + / 2 r 


2 .... Only lines from n u < 10 with 1.05 > /(r 388 9) > 0.95 



over the range T3889 = [0, 100] are included. 
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Table 4. Limited Range Fitting Formulae For Helium Abundance Including Optical 

Depth Effects a 



Line 

(A) 


A 


B 




a 


bo 




6l 


b 2 


Max Error 

(%) 


3889 


0.904 


-0.173 


-5.4 x 10~ 4 


-1.06 x 10" 1 


5.14 x 10" 


-5 


-4.20 x 10~ 7 


1.97 x lO" 10 


0.8% 


4026 


4.297 


0.090 


-6.3 x 10~ 6 


1.43 x 10~ 3 


4.05 x 10" 


-4 


3.63 x 10~ 8 




0.2% 


4387 


16.255 


0.109 


-3.7 x 10~ 6 












0.2% 


4471 


2.010 


0.127 


-4.1 x 10~ 4 


2.74 x 10~ 3 


8.81 x 10 


-4 


-1.21 x 10~ 6 




1.0% 


4686 


0.0816 


0.145 














0.3% 


4714 


19.555 


-0.461 


-13.5 x 10~ 4 


5.77 x 10~ 2 


-7.12 x 10" 


-3 


-3.10 x 10~ 5 


3.99 x 10~ 8 


1.9% 


4922 


7.523 


0.150 


-1.9 x 10~ 4 












0.6% 


5876 


0.735 


0.230 


-6.3 x 10~ 4 


4.70 x 10~ 3 


2.23 x 10 


-3 


-2.51 x 10~ 6 




1.2% 


6678 


2.580 


0.249 


-2.0 x 10~ 4 












0.3% 


7065 


See text 


















7281 


13.835 


-0.307 


-7.0 x 10~ 4 












2.0% 



a Helium number abundance is given by y = (In ne / '1 ' Hp)F(n e ,T) / f{n e ,T,T3gg<j), where Iu ne /I[{f3 is the observed line ratio, 
F(n e ,T) = At B ° +Bl ™" (from Paper I), and the optical depth correction factor is given by /(n e , T, T3889) = 1 + (r3889/2)[a + 
(60 + bin e + fe2«e)i]- These functions apply only for cases in the range n e = 1 — 300 cm" " 3 and T = 1.2 — 2.0 X 10 4 K. The 
maximum uncertainty over this range is noted. 



